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Abstract. We find convergent double series expansions for Legendre's third 
incomplete elliptic integral valid in overlapping subdomains of the unit square. 
Truncated expansions provide asymptotic approximations in the neighbour- 
hood of the logarithmic singularity (1,1) if one of the variables approaches this 
point faster than the other. Each approximation is accompanied by an error 
bound. For a curve with an arbitrary slope at (1, 1) our expansions can be 
rearranged into asymptotic expansions depending on a point on the curve. For 
reader's convenience we give some numeric examples and explicit expressions 
for low-order approximations. 
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1. Introduction. Legendre's incomplete elliptic integral (EI) of the third kind is denned by 
formula 13.6(3)] 

A 

U(X,V,k)= [ ; ^ ■ (1) 

J (l + vt 2 )^/(l-t 2 )(l-k 2 t 2 ) 

It is one of the three canonical forms given by Legendre in terms of which all elliptic integrals can be 
expressed. We will only consider the most important case 0<fc<l,0<A<l and v > — 1. Two 
other canonical forms can be expressed in terms of n(A, is, k): setting v = gives the first incomplete 
elliptic integral, while setting v = — k 2 gives the second plus an elementary function (see formulas 
13.6(11)-13.6(13)]). When A = 1 the rhs of Q is called the third complete elliptic integral. 

Series expansions, asymptotic approximations and inequalities for the third incomplete elliptic 
integral have been studied by many authors. Major contributions were made by Radon [Sj, Carlson 
PJ|2], Carlson and Gustafson 3 and Lopez jUlE]. Carlson showed in that H~(A, u, k) can be expressed 
in terms of Lauricella hyper geometric function Fp of three variables (see (I]). In the same paper he 
noted that one can derive rapidly convergent expansions by first expressing Legendre's incomplete Els 
in a different form. This form later became known as symmetric standard Els. Today, most authors 
consider asymptotic approximations for these symmetric Els. Since Legedre's incomplete elliptic 
integrals are connected with symmetric Els by certain simple relations (see E] ) , it is possible to 
reformulate expansions presented here in terms of symmetric elliptic integrals. 

In this paper we give convergent double series expansion for the third incomplete elliptic integral 
valid in two overlapping subregions of the unit square [0, 1] x [0, 1] in (k 2 , A 2 )-plane defined by 

Region I = {k 2 , A 2 G [0, 1] : A 2 < 1/(2 - k 2 ) & (1 - A: 2 )A 2 /(1 - A 2 ) < 1}, (2) 



Region II = {k 2 , A 2 G [0, 1] : A 2 > 2 — 1/A; 2 (1 - A 2 )A: 2 /(1 - k 2 ) < 1}. (3) 
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Combined these subregions cover the unit square completely. Truncating our expansions one gets 
asymptotic approximations for (1 — A)/(l — k) — > or (1 — A)/(l — k) — ► oo, i.e. for curves with 
endpoint on the top or right side of the unit square or with endpoint A = k = 1 and horizontal or 
vertical tangent at that point. In all cases we give explicit bounds for the remainder term. 

For a given smooth curve approaching the singular point (1,1) from any direction our double series 
expansions can be rearranged into asymptotic expansions, where a point on the curve plays the role 
of an asymptotic parameter. In the last section of the paper we illustrate our results by numerical 
examples. Rather good precision (with less than one percent relative error) is reached in most cases 
already by the first asymptotic term. For reader's convenience we also give explicit expressions for 
first several terms of each expansion. 



2. Partial fraction decompositions. For a non-negative integer n and a = 0, 1 define 

6 a (t) = (I -fr™ + (-l) a (l+t)~ n w a (t) = l K 1- *) n -(-l) a (l+0 n -2a]/n, n > 0, , . 
W) U V +{ i) U + tj , Wtj | (-l)°ln(l + t)-ln(l-t), n = 0. 1 j 

The second definition is consistent in the sense that lim ip%(t) = (—1)° hi(l + 1) — ln(l — t), a = 0, 1. 
Our first lemma is straightforward. 

Lemma 1 The functions and are related by: 

A 

J </>%(t)dt = ^_ 1 (\), n = 1,2,..., a = 0,1. (5) 
o 

Lemma 2 Let n < 2j be non-negative integers. Then the following partial fraction decomposition 
holds true: 

j r ^~ = &[imt) + (?[2\m) +■■■+ ®mm)> ( 6 ) 

where a = nmod2, P^\j\ = 2~' J and the numbers /3™[j — k], k = 1, 2, . . . ,j — 1, are found from the 
following recurrence (here (£) = when k > n): 

~ k] = ^ Q " E(-2r fe ( fc { -) k = 1, 2, . . . ,j - 1, (7) 

or explicitly from 



OT-fc]=2-^E(-2) l 0( fc + ]_\ 0<k<j. (8) 

Proof. For computing the coefficients of the partial fraction decomposition we use the trick 
routinely applied for calculating residues. Writing 

r _ [l] ff [2] fffl 

+ l-t (i-i) 2 (i-t) j (i + J '' 

where 

+ 77^ TT7 H h 



(i + ty i + t (i + t) 2 {i+ty : 

and multiplying throughout by (1 — t 2 )- 7 , we get: 

t n = (i + ^(^[j] + ^ - i](i - 1) + • • • + otki - ty- 1 ) + (i - tyR(t). 

Setting t = 1 immediately reveals = 2~ J . Differentiate fc times and set t = 1: 

HfSi = E (*) [(i + ^'iftf^b- - *K-iM fc = i, 2, . . . , j - 1. 

i=0 ^ ' 



We have: 



[tJ l*=i"\0, fc>n ' Ki + tj J l*=i " (j-jfe + t)! 



2 J 



-fc+i 



Thus the numbers (3p[j — fc] can be found via the following recurrence (starting with ffplj] = 2 3 and 
counting downwards): 



^ " *'<-«***' " " g (*) (P^^-^b - SK-DM. * = 1,2, 1. 

A simple rearrangement of this formula yields which can also be written as 

- / a \ f 



I>^>-^(I) 



This recurrence can be solved for f3™[j — i] using the inversion formula (see [2]) 



,n — k I ^— ' V n — k 

k=o K 7 k=o v 

which leads to ©. □ 

In the sequel we will only need even n therefore we put <p n (t) = 0° (t) and ip n {t) = ipn(t) to simplify 
notation. 

Lemma 3 The following partial fraction decomposition holds true for j = 0, 1, . . .: 

t 2j _ f (^)Vj(0 f Y 

(l + „t 2 )(l-t 2 y+ 1 (l + z/)i+l(l + ^ 2 ) ^ (l + z/)^ 1 ^ + J ' l j 

where the numbers a n)in are expressed in terms of the numbers (3p[i] defined by ((7J) or (0) &y 

«n, m = (-i) m 4+r m - 4 N, n = i,...,j+i ; ™={J;;;; „>i ; ai -° = 5- (10) 

Remark 1. This lemma contains two statements: one is formula (jlUjl for the coefficients of expansion 
and the other is that the numbers a n ^ m are independent of j. That is when j increases we do 
not need to update all the numbers a n m , instead we keep all previously calculated numbers and 
complement them with new j + 1 numbers with indices summing up to j + 1: n + m = j + 1. 
Proof. We will use induction in j. For j = we can verify directly that: 



1 v 1 

+ 



(l + ^2)(1_f2) (l + I/ )(l + 1/t 2) 2(1 + 1/) 

Suppose now that © holds for a fixed j. Then for j + 1, 

/-'•' ' 2 



1 1 

+ 



1-t 1+t 



{I + vt 2 ){l - t 2 )i+ 2 {I + - t 2 )i+ 2 {l + u){l + ut 2 ){l-t 2 )i+ 1 ' 

Now using © with n = 2j, m = j + 2 for the first term and © for the second gives after collecting 
coefficients at (j>i(t): 

t 2 i+ 2 {-\y+ l u 



!¥_mn -i- u ) j+1 



(1 + ut 2 )(l - t 2 )i+ 2 (1 + z/)i+ 2 (l + z/i 2 ) + 
(-l)i+i( aij0 + a 1)1 (l + i/) + --- + a 1)J -(l + z/)i)+^ 2 [l](l + 

+ (TT^ 0l( } + 

(_l)i ( a2>0 + a21 (i + „) + ... + a 2j _!(l + i/)^ 1 ) + $J 2 [2](1 + r/)^ 1 

+ (TT^V Mt) + 



+ 



(TTW 2 * i+l(t) + i + * 0J+2(t) - 



Now define 



a 1|i+1 = (-1)'' +1 $J 2 [1], a 2tj = (~iy^ 2 [2], a j+2>0 = P*i 2 \j + 2). 



?2i 



This shows that expansion for j + 1 has the form Q with numbers a n ^ m given by ()1U() . □ 
Formula © combined with Lemma pleads to the following evaluation 



A 

I 



0dt 



/J} 

i + ^y+i arctan ( A ^) + ^d jtn+ i{v)iij n (\) 



(l + vt 2 ){l ~t 2 y+ X + 



(11) 



n=0 



where we introduced the notation 



dj in+ i(v) 



(-1) 



j-n 



(1 + 



^a n+ i,i(l + f) 4 , n = 0, 1,2, . . . , j. 



(12) 



Here and henceforth y / i/arctan(Ay / z/) will be understood as (V — ^)arctanh(A\/ — ^) for negative v. 
Another way of decomposing the left hand side of Q is given in the next lemma. 



Lemma 4 The following expansions hold true for j = 0, 1, . . .: 







{I + ut 2 ){l - t 2 )i 



+i 



1+2/ 



1 + ,-/ 2 + S (l + // )./-.'H-2 (| _ ,2,, 
n=l ' 



and 



(-i/uy 

l + vt 2 ~ 1 + ut 2 



E(-v 



(13) 



(14) 



i=0 



Proof. Both expansions can be verified by using the summation formula for a finite geometric 
progression. □ 

Euler's integral representation for the Gauss hyper geometric function 2 Fi (after simple variable 
change) combined with ® and (JHJ) gives 







0dt 



(1 - t 2 ) n 2j + 



- 2 Fi(n, j + 1/2; j + 3/2; A 2 ) = £ /^^(A). 



(15) 



Lemma |1] and formula (|15[) provide two alternative ways to evaluate the integral in (|11|) . We have 
A „. ■ , -, A 



t 2j dt 



o 



{I + vt 2 )(l - t 2 )i+ l 



l + v 



j+1 f t 2 Ht X 2j+1 ^ ^-»+ 1 2 F 1 (n,j + l/2;i+3/2;A 2 ) 

+ ^-^2^ n i .AT-r,.4-2 ( 16 ) 







1 + z4 2 2j + 



n=l 



(1 + i/)i-«+ 2 



f_lV fu (— IV" J_1 A 2i+1 J+1 jj-n+l n 
= T ^XT arctan(A^) + 77^ ^rV(-i/) i+1 - +V 7 — — V /^[#i_i(A). (17) 

j=0 n=l i=l 

3. Expansions for the third incomplete elliptic integral. The following relation will be of 
great help: 

II(A, v, k) = U(u, k) - ^===11 (-v/(l + u), Vl-A 2 , ^-k 2 /{l-k 2 )) , (18) 

(1 + u)\J\ — k z \ / 

where n(i/, /c) is the complete elliptic integral of the third kind. Here we choose the branch of the first 
square root that is positive for positive values of 1 — A 2 . The choice of the branch of the second square 
root is immaterial since n depends on the squared second argument only. This relation can be easily 
verified by representing the integral over (0, A) from Q as the difference of integrals over (0, 1) and 
(A, 1) and introducing the new integration variable u 2 = 1 — t 2 . 



Theorem 1 Let k,X 6 Region I and v > —1. JTien /or any 'positive integer N the following 
expansions hold true 



N-l 



j=0 



(-i) J V^ 

(1 + 



arctan(A^) + ^ dj >n+1 (v)i() n (A) 



n=0 



+i? liA r(A;,A,i/), 
(19) 



and 



H(A, „,fc) = £ - k 2 y { arctan(A^) 

j! [(l + i/)^ 



+- 



A 



2i+l j±i 



2? + I + E n + y )i-n +2 E M^iW f+«i,*(M,*'), (20) 

n=l 



v 7 i=0 n=\ v ' i=l 

where the functions ^n(A) and dj )n (v) are defined by Q) and (|12|) . respectively, and the remainder 
term satisfies 



\Riji{k,\,v)\ < 



(1/2)jvA 



2min(l,l + i/)iViV! 



(l-fc 2 )A 2 
1-A 2 



JV 



(21) 



Proof. Put k' 2 = 1 — k z and calculate using binomial expansion and termwise integration: 



n(A, v, k) 



dt 



(1 + vt 2 )y/(l -t 2 ){l-k 2 t 2 ) 7 (1 + ^ 2 )(1 - f 2 ) I 1-t 2 



dt 



1 + 



*/vv 1/2 



dt 



' N-l 



{1 + vt 2 ){l - t 2 ) 

N-l 

= E(-d 



v (1/2), « 2 H 2 i " - (1/2), ^gi 



3=0 



t 2 -?'dt 



(l + ^ 2 )(l-t 2 )i+ 1 



+ i?i ) jv(/s,A 3 i/) ) 



where the remainder term is given by 

^^^,^ = ^(-1/^(1-^7- 



t 2j 'dt 

(1 + z/t 2 )(l - i 2 )J +1 ' 



(22) 



(23) 



Using representations (fTTj) and (|T7|) for the integral in (|22j) we get (|T9*jl and (|2()|). respectively. To 
prove the error bound (|21j) we note that the series (|23j) has alternating signs. The following estimate 
shows that the terms in (j23j) monotonically decrease in absolute value: 



(i + D! ( } 



t 2 i+ 2 dt 



(l/2),-(l/2 +j) /■ (l-A; 2 )^ * 2 (l-fc 2 ) 



(l + z4 2 )(l-t 2 )i+ 2 + 1) 7 (l + ^ 2 )(l-t 2 )i+! (l-t 2 ) 



dt 



< 



(1/2), 



(1 - A: 2 )' 



f 2; 



x\i-k 2 ) dt<m± 



(l + vt 2 )(l-t 2 )i+ l (1-A 2 ) 



-(l-fc 2 y 



The last inequality is due to (J2J). Hence we are in the position to apply the Leibnitz convergence test 
and the remainder term does not exceed the first term in (|23[): 



|fli^(M,*)l<^(l-*Y 



t 2N dt 



< 



(l/2) N (l-k 2 ) 



N 



(l + vt 2 )(l-t 2 ) N + 1 ~ jV!min(l,l + i/) J (l-t 2 ) 



t 2N dt 



\N+1 



The integral on the right hand side satisfies the following asymptotically precise (as A — * 1) estimate 



r t 2N dt X 2N+1 

= J (j—pym ^ 2N(1 - A 2 )^ = /2 ^' ( 24 ) 
o 

which is valid for all A G (0, 1) and N > (not necessarily integer). Indeed, /i(0) = /2(0) = and 

2N , f , 

tAy = 2JV + 1 _ A2 < 1 - WD- 

This estimate immediately leads to (|21j) . □ 

Remark 2. The error bound (|21|) shows that expansions ()19|) and (|20j) are asymptotic for A — ► 1 
with either constant or —* 1 in a way that (1 — A)/(l — &) — > 0, so that we can approach the 
singular point (1, 1) along a curve in (k, A) plane having infinite slope at (1, 1). These expansions can 
be transformed into asymptotic expansion for k, A — ► 1 along an arbitrary curve if an ordering can 
be introduced into the matrix of functions forming the expansion that converts this matrix into an 
asymptotic scale. Indeed, we have the matrix 

V>o(A) 
I 

arctan(A- v /Z') 

I 

(1 - k 2 ) arctan(A v / ^) (1 - k 2 )^{\) <- (1 - k 2 )^x(X) 
(l-£: 2 ) 2 arctan(A^) (1 - k 2 ) 2 ^{\) «- (1 - £; 2 ) 2 Vi(A) <- (1 - fc 2 ) 2 V 2 (A) 

■!■ "1* *!■ 



Arrows indicate the ordering that is true no matter from which direction we approach the point (1, 1). 
For instance, (1 — k 2 ) 2 ipi(X) = o((l — /c 2 )V>i(A)) as k, A — > 1. However, unless we specify a curve in the 
(k, A) plane this is only a partial ordering. For instance, we do not know how to compare (1 — /c 2 ) 2 V'i(A) 
and (1 — A; 2 ) 3 V'2(A). Choosing a curve possessing a tangent at the point (1, 1) introduces the complete 
ordering into the above matrix. Rearranging terms in the expansions (|19j) or (|20j) according to this 
ordering converts it into the asymptotic expansion for the specified curve. If the slope of our curve 
at (1, 1) is strictly between and vr/2, then all the diagonals of the above matrix will be of the same 
asymptotic order and the order will be decreasing as we move away from the main diagonal. 
Remark 3. The error bound (|21|) shows that expansions (j!9ft and (|2()jl are convergent for k, A in 
Region I. Taking N = oo in (|19|) and (|2U|) and summing up we obtain convergent double series 
representations 



n(A, v, k) = J^{» + k 2 )- 1 ' 2 arctan(Av^) + £ £ f ^ i.n+l Wn(A), (25) 



n<A ' "' k) = J (TT^VTT^ + g S RWWTTKTTT^^ 1 (», ; + 1/2U + 3/2; A ), 


(26) 

where = u(l — k 2 )/{l + u). The first series requires additional condition —k 2 < v for convergence. 

To derive an expansion valid in Region II we start with an expansion valid in a neighbourhood of 
(0,0) given by Carlson pQ: 

00 \2m+l m (-1 ln\ 

H(A, v, k) = £ A-^ ^(-^-li-A^-n, 1/2; 1/2 - n; k 2 ). (27) 

m=0 n=0 



Using 



(1/2) 



— 2 F l {-n, 1/2; 1/2 - n; k 2 ) = 2 F x {-n, 1/2; 1; 1 - k 2 ) 



(28) 



which is a limiting case of the well-known analytic extension formula for 2 Pi (see jlj formula 2.10(1)]), 
we obtain the expansion 



^2m+l 



n(A, u,k) = J2 ^—pi Y,(-^ m ~ n * F ^- n > V2; i; 1 - fc 2 



(29) 



m=0 



n=0 



which is valid for \v\\ 2 < 1 and \k\ < 1/A as shown by Carlson. An application of (|TH|) leads to 
Theorem 2 For /c, A S Region II, ^ in the range 

((l-A a )M)/(* + l)<l, (30) 
and positive integer N the following expansion holds true 

2\m _™ ^-« 2j p i (_ n) 1/2; 1; (1 _ k 2 )- 1 ^ 



U(X,u,k) = U(u, k) 



1-A 2 '^ (1-A 2 ) ^ 

m=0 n=0 



(1 + v)m-n+l 



+ R 2jN (k,v, A), 
(31) 



wi/i i/ie error bound given by 

\R 2 ,N(k,v,\)\ < 
where f(\,k,v) is defined by 



(1 _ \2)N+l/2 

2N + 1 



max 



A- 2 



1 - A: 2 ' v + 1 ' 



-1 JV+1/2 



/(A,*!,!/) 



(32) 



/(A, k, v) 



{l+u)k 





)IH 


1 




V(i-fc 2 )(i+v 


)IH 


1 


1 - 


(i+^)Vi-fe 2 


1 





(1-A 2 )fc 2 
1-fc 2 

1 - 
1 - 



-1 r 



1 o+i)fc 2 



(1-A 2 )H 



V(i-fc 2 )(i+^)ki 







(1 


-A 2 )k| 




l+v 


A 2 ) 


HI 


-1 








(1 


-A 2 )M 



-1 



-1 



max 



k 2 IH 



1 



(H-l)fc 2 
IH(i-fe 2 ) 



-1 



k 2 



> l; 



IH ^ k 2 

u+l T-k? 



1 l+v 
1 IH 



1 _ _M_ 

i+v 



(i-a 2 )H 



H_ >> 1 > k 2 . 



+ iV), max( T ^,lj=^ T . 



IH 



(33) 



Proof. To derive the form of (|31|) expand the second term on the righthand side of (|18|) into series 
(|29jl . Carlson's conditions |^|A 2 < 1 and < 1/A written for the parameters of this second term 
become (0) and (|30|) . Now consider the remainder term 



R N (k, u, A) 



1 



l + v 

(1 - A 2 ) m 



1-A 2 



1-k 2 



1/2 



R N (k,v,\) = Y, fclfE ; 



2m + 1 ^ \ 1 + 

m=N n=0 



R N (k,f, A), 

2 Fi(-n,l/2;l;(l-A ; 2 r 1 ). 



Change the order of summation to get 



»/, A) = J] J] ^ ; !„.. AnN ; ; ( T~7"77 ) 2^l(-n, 1/2; 1; (1 - fc 2 ^ 1 

n=0 s=0 ^ 



(s + max(n,iV)) + 1 Vl + ^ 



Dl-A^VlUlAlia-fc 2 )" 1 ) 



n=0 



, I - A>] £ (1-A 2 )>/(1 + ,))- (34) 



l + i/ 



s=0 



2(s + max(n,iV)) + 1' 



Since (L > -1/2) 



we may write 



E 



x 



L-l/2 r fL-l/2 



2s + 2L + 1 



1 -t 



R N (k,u,X) = - 



1 + u 



(1-A> 



1/2 



v^(l-t) 



E 

n=0 



l + V 



t ra& < n ' N \F 1 (-n, 1/2; 1; (1 - A; 2 ) -1 ). 



The series in (|35j) can be broken into two parts: 
JV-l 



n=0 



2 F 1 (-n,l/2;l;(l-A; 2 )- 1 )+ ^ 

n=iV 

To give a bound for we will need two ingredients. First is the estimate 
la^iC-n, 1/2,1,(1-*) )|< <fc 2 <l/2, 



(35) 



f^-n, 1/2; 1; (1-fc 2 )- 1 ). (36) 



(37) 



which is obtained by writing 2F1 as the Legendre polynomial |10| formula(7.3.176)]: 



2 F 1 (-n,l/2;l;(l-k 2 )- 1 ) 



r —k 2 - 


f / 


Pn 1 


l-k 2 _ 



1-2A; 2 



2a/-/c 2 (1 - A; 2 ) 



A' 2 



1-A; 2 



P 



-i(l-2k 2 



(38) 

and applying the inequality (valid for all complex z) \P n {z)\ < \z + y/ z 2 — l| n , where the branch of 
the square root is chosen so that \z + \J z 2 — 1| > 1. This inequality follows immediately from the first 
Laplace integral for P n {z) (see jU formula 3.7(6)]). The second ingredient is the asymptotically precise 
(as x — ► 0) estimate (verified by comparing derivatives of both sides and noticing that /(0) = g(0)): 



/(*) 



t M dt 



,M+1 



< 



1-t ~ (M + l)(l-x) 



(39) 



valid for M > — 1 and < x < 1. To prove (|32j) - (|33j) we need to consider the following cases: 
a) 1/2 < A: 2 < 1, (1 + v)k 2 /{\v\{l - k 2 )) > 1; b) 1/2 < A; 2 < 1, (1 + v)k 2 /{\v\ (1 - A; 2 )) < 1; 
c) 1/2 < A; 2 < 1, (1 + v)k 2 /(\v\{l - k 2 )) = 1; d) < A: 2 < 1/2, v > -1/2; 
e) < A; 2 < 1/2, -I < v < -1/2; and f) < A; 2 < 1/2, v = -1/2. 

Take a). Inequality (|37j) gives for the first sum in (|3*d]> : 



N-l 

E 

ra=0 



1 + V 



2 F 1 (-n,l/2;l;(l-fc 2 )- 1 )| < 



t(l + i/)A; 2 



— A; 2 ) 



v 



[(»/ + l)^]/[H(l -*»)]-!• 



Similarly, the second sum in (|36|) satisfies 



E 

n=N 



l + V 



t n 2 F 1 (-n,l/2;l;(l-k 



2\-l- 



< 





t(l 


+^)fc 2 


AT 




v(l-fc 2 ) 




1 - 


t(l+z/)fc 2 


;/(l-fc 2 ) 



< 



t(l+v)k 2 



N 



-. _ (1-A 2 )fc 2 ' 
(1-fc 2 ) 



The last inequality is due to the fact that t £ [0, (1 — X 2 )\u\/(u + 1)]. Hence the complete sum (|3*6*)l is 
estimated by 



E 

n=0 



l + V 



t max{ - n ' N+x \F x {-n,l/2; 1; (1 - A: 2 )" 1 ; 



< 



t(l + v)k 2 



v{l-k 2 ) 



N 



fc 2 (l-A 2 ) 
(1-fc 2 ) 



u+1 _ 1 
IH(1-A 2 ) 1 



(i/+l)fc 2 _ -, 



(1-A 2 )fc 2 
1-fc 2 



Substituting this estimate into and using (|3"H|) we get: 

l + i/ 



\R N (k,v,X)\ < 



(2N + l)\u\ 



(l-A 2 )fc 2 nA 
(1-A: 2 ) _ 



k 2 /{l-k 2 



1 



1 



(1-A 2 )fc 2 



which is equivalent to (|32[) plus first line of ((SSI) ■ Proofs in all other cases require only minor modifi- 
cations. □ 

Remark 4. Error bound (|H2|) shows that expansion (j31j) is convergent and asymptotic as — > 1 with 
either constant A or A — > 1 in a way that (1 — k)/(l — A) — > 0, so that the point (1, 1) is approached 
along a curve in (k, A) plane having zero slope at (1,1). It can be rearranged into an asymptotic 
expansion for k, A — > 1 along an arbitrary curve as explained in Remark |2J 

Remark 5. Expansion 1)31 J) can be derived without any use of Carlson's expansion (|27|) by expanding 
the rightmost term in ()18|) into binomial series, similarly to the proof of Theorem^ 
Remark 6. The complete elliptic integral n(^, k) in (j31j) can be approximated using the following 
asymptotic expansion valid for k — > 1: 

n(^ fc)= £ (1/2) ^ 1 / 2)n (l - k 2 ) n m + n) - ^(1/2 + n)-\ ln(l - fc 2 )) 2 *i(l, 1/2 + n; 1/2; 



n=0 
M-l 

E 

n=0 



(l/2) w (l/2) n 
2fn!) 2 



(1 - fc 2 )^ ' 1 ' ' ^!, 1/2 + n ; 1/2; -i/) + ((1 - k) M In ^-L^ 



(40) 



Here 

2^(0,1,0,0) is 

the derivative of 2F1 m the second parameter. This expansion can be obtained by 
writing II(i/, k) as Appell's hypergeometric function F\ (see pQ), expressing Fx as a series in terms of 
Gauss hypergeometric functions 2-F1 and applying [4, formula 2.10(12)] to 2-^1- 



4. Results of computations. In this section we give several examples of computations with the 
expansions obtained above. Denote by Iljv(A, u, k) the iV-th order approximation in i.e. the right 
hand side of ()19|) without the remainder term R\^{X, u, k). Then: 



Ui(X,u,k) 



1 



2(l + i/) 1-A 



1 + A v/z/ arctan (A-i/z/) 
In r H — 



l + i/ 



U 2 (X,u,k) = IL 1 (X,v,k) 



1 ~ k L /-, nf- 1 , i + A 

; -7 -r U + (1 - A In - 

l + i/ 1-A \ v >v + l 1-A 



1 + 1/ 



arctan (Ay^) 



1-A 

TTa 



2\2 



n 3 (A,i/,fc) =n 2 (A,z/,fc) + 



(i-fe 2 ) 



128(l + i/) (1-A) 2 



1 - 2A 



l + i/ 



+ 1 



1 - A 

TTa 



in 



l + A 



1- A + 



1 (1-A) 2 + 16 



(1 — A) 2 ^fv arctan (X^/v) (1 — XY 



:i+^ 



a+Ar 



A 


fe 


1 - k 
1 - A 


U(v,X,k) 


niO/,A,fc) 


relative 
error 


relative 
error 
bound 


U 3 (X,u, k) 


relative 
error 


relative 
error 
bound 


.5 


.9 


.2 


.37138 


.37409 


-.730xl0" 2 


.213X10- 1 


.37138 


-.634xl0~ b 


.178x10-* 


.6 


.99 


.025 


.41973 


.42022 


-.117xl0~ 2 


.400xl0~ 2 


.41973 


-.306xl0" y 


.104xlO- H 


.75 


.999 


.004 


.48662 


.48673 


-.218xl0~ ;i 


.990xl0~ ;i 


.48662 


-.289xl0" y 


.136xl0~ 8 


.9 


.99999 


.0001 


.57202 


.57203 


-.553x10^ 


.335xl0~ 4 


.57202 


-.784xl0~ 14 


.508xl0~ 13 



Table 4.1. Numerical examples for the approximation (|19|) with k — > 1, X — ► 1, (1 — k)/(l — X) ^ 
and v = 7. The sixth and the ninth columns represents the relative errors i?i i Ar(A, z/, k)/U(X, u, k) in 



(j!9j) . The seventh and the tenths columns represent relative error bounds as given by the rhs of (|21[) 
divided by n(A, u, k). 

For k and A satisfying Q denote by 11^ (A, v, k) the iV-th order approximation from 1)31(1 . i.e. the 
right hand side of 1)31(1 without remainder term i?2,7v(A, v, k). Then: 



IIi (A, i/, A;) = U(u,k) 



1 



1-A 2 \ 1/2 
1 +/' \ l-k 2 ) ' 



_ 1 f\ _ \2\ 3 / 2 ( i _ u2 , 

u 2 (x,u,k) = u 1 (\, v ,k)--i^— ¥ j [j^ln- 



i i 



1 + 1/7 2 l + i/ J ' 



n 3 (A,z/,/c) 



n 2 (A,i/,fc)-i (i- 



A 2x5/2 



fc 2 J 



(1-fc 2 ) 2 
1 + 1/ 



1 + 



+ 



1 + 1/ 
1 1 

~2 l + i/ 



(l + i/) 2 
A; 2 



2 + 



l + i/ 



3 

+ :r 



A 


k 


1 - A 

1 - k 


n(A,i/, &) 


fii(A,i/, k) 


relative 
error 


relative 
error 
bound 


n 3 (A,!/, k) 


relative 
error 


relative 
error 
bound 


.9 


.5 


.2 


.50760 


.51315 


109x10^ 


.867X10- 1 


.50770 


-.212xl0" 3 


.135x10^ 


.99 


.6 


.025 


.56514 


.56530 


-.287xlO~ 3 


.238xl0~ 2 


.56514 


-.606xl0~ 7 


.403xl0- fe ' 


.999 


.75 


.004 


.60555 


.60556 


-.682xl0~ b 


.376xl0- 4 


.60555 


-.138xl0- lu 


.106xl0" y 


.9999 


.8 


.0005 


.62453 


.62453 


-.153xl0~ b 


.110xl0" b 


.62453 


-.888xl0" lb 


.597xl0~ i;5 


.999999 


.95 


.00002 


.71429 


.71429 


.172xl0~ 8 


.540xl0~ 8 


.71429 


.190xl0~ 18 


.793xlO- iS 



Table 4.2. Numerical examples for the approximation l|31|) with A — > 1, k — ► 1, (1 — A)/(l — k) — > 
and v = 7. The sixth and the ninth columns represents the relative errors R2 jv(A, v i &)/n(A, u, k) 
in ()31(), The seventh and the tenths columns represent relative error bounds as given by the rhs of 
divided byU(X,u,k). 
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